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ABSTRACT 



The vicinity of the supermassive black hole associated with the compact radio 
source Sagittarius (Sgr) A* is believed to dominate the observed emission at 
wavelengths near and shorter than ~ 1 millimeter. We show that a general 
relativistic accretion flow, heated via the plasma wave heating mechanism, is 
consistent with the polarization and recent mm- VLBI observations of Sgr A* for 
an inclination angle of ~ 45°, position angle of ~ 140°, and spin < 0.9. Structure 
in visibilities produced by the black hole shadow can potentially be observed 
by 1.3 mm- VLBI on the existing Hawaii-CARMA and Hawaii-SMT baselines. 
We also consider eight additional potential mm- VLBI stations, including sites 
in Chile and New Zealand, finding that with these the basic geometry of the 
emission region can be reliably estimated. 

Subject headings: Galaxy: center — black hole physics — accretion: accretion 
disks — sub-millimeter — techniques: interferometric 



1. Introduction 



The compact radio source Sagittarius (Sgr) A*, residi ng at the Galactic center, is be- 



lieved to be the best candidate of suppermassive black hole (jSchodel et al.ll2002l : iGhez et al. 



20051 ). Based on observations of the orbiting stars, the mass of its central object is measured 
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as over 4 million M & f lGhez et al.ll2008l ; iGillessen et al.l 120091 ) . However, it remains to be 



conclusively proved that Sgr A* is a black hole, and the explanations of the broad-band 
emission are even controversial. Various morphology of the emission region, torus, jet, and 
their combinations, are provided by different models. The observations at millimeter or 
longer wavelengths cannot distinguish them due to the interstellar scattering. However, the 
interstellar scattering is becoming less severe as the observational wavelength decreases to 



(Bower et al. 


2004; 


Shen et al. 


2005) 



trinsic structure of Sgr A* would be dominant, i.e., the sub-millimeter VLBI promises to 
produce unadulterated images of its intrinsic structure. However, lack of baselines results in 
poor uv-coverage, making the imaging impossible. Thus, theoretical modelling of Sgr A* will 



be cr itical to both planning and interpreting sub-millimeter VLBI observaions (IHuang et al. 
2007h . 



Recently, iDoeleman et al.l (120081 ) reported observations at 1.3 mm with a very- long- 
baseline interferometry (VLBI) array consisting of the Arizona Radio Observatory Sub- 
millimeter Telescope (SMT), Combined Array for Research in Millimeter- wave Astronomy 
(CARMA), and James Clerk Maxwell Telescope (JCMT). Robust detections are obtained on 
baseline CARMA-SMT at ~ 500 km and baseline JCMT-SMT at ~ 4500 km. An upper limit 
is yielded on baseline CARMA-JCMT at ~ 3000 km. Since these visibility measurements 
are too few to produce an image, a circular Gaussian modelling results in a full width at half 
maxi mum (FWHM) of ~ 37tio/uas, with the scattering effects removed (IDoeleman et al. 



20081 ). However, also pointed out by them, this size is somewhat unexpected, because it 
is smaller than the 'black hole shadow' size of ~ 50/ias regardless of the black hole spin 



( jFalcke et al.l 120001 ) . It might be possible that the emission really comes from the inside of 
black hole shadow. On the other hand, a structure with a central black hole shadow itself, 
deviating from the Gaussian, might be implied by the visibility measurements. 

In this paper, we present a general relativistic accretion flow model with plasma wave 
heating mechanism (SecJS]). The polarization observations can be well explained with ap- 
propriate viewing angle and other model parameters set for arbitrary black hole spin. We 
then tend to use eight potential stations for future 1.3mm VLBI measurements to perform 
visibility analysis (SecJ3|). We perform a visibility analysis using eight potential stations for 
future 1.3 mm VLBI measurements. We test the images with recent 1.3 mm VLBI measure- 
ments and investigate the implications for the properties of the black hole and accetion flow 
(SecHJ). Comparisons to the earlier work and evaluations on the potential (u, v) coverage 
are given in the final Section. 



- 3- 



2. General Relativistic Keplerian Accretion Flow 



We follow the m agneto-rotational-instability (MRI) driven accretion flow (IMelia et al. 



200ll ; iLiu et al.ll2007l ) where the primary mechanism for generation of turbulence and vis- 
cous stress in accretion flow s supported by magnetohydrodynamics (MHD) simulations 
( IBalbus fc Hawleylll99ll . Il998l ). When deriving the basic equations of the general relativistic 
accretion flow, we use the energy momentum tensor T^ v given as T^ w = n(m p + m e )r)u tl u' / + 
pgV _|_ fiw _|_ qV u v _|_ wnere n ^ t^", p ; 7^ mp anc l me are the rest number density, 

four- velocity, heat-flux four- vector, viscous tensor, pressure, relativistic enthalpy, the mass of 
the proton and the mass of the electron, respectively. The basic equations for the relativis- 



tic hydrodynamics are then derived from the baryon-mass conservation V, 



[nw 



0, the 



energy-moment um conservation y v T^ v = and an equation of state. Detailed derivations 
can be found in iTakahashil (120071 ). In this study, the Boyer-Lindquist coordinate is used for 
the description of the Kerr metric. 



Following ILiu et al.l (120071 ) and iHuang et al.l (120081 ). we introduce model parameters as 



the ratio of the total stress to the magnetic field energy density f3 v , the ratio of the magnetic 
field energy density to the gas pressure j3 p , dimensionless constant for electron heating rate 
C\, constant mass accretion rate M, inclination angle i, and position angle B. An additional 



meas ured from the orbit of the short-period star SO-2 for the distance of 8.0 kpc (IGhez et al. 



1994 ; lAbramowicz et al.lll996l ) 



(Ghez et al. 


t\ ( 


Lasota 



—n(m p + m e 



^3/2 A l/2 7 3 



dVt j dr, 



where A = r 4 + r 2 a 2 + 2Mra 2 } A = r 2 — 2Mr + a 2 , 7 is the Lorentz factor, = 1$ ju 1 = d(j)/dt 
is the angular velocity of the accretion flow, and v is the kinematic viscosity coefficient. We 
then assume the a-viscosity as 



V 



n 



m p + m e )(— rdfl/dr) ' 



(2) 



where the product of (3 V and j3 p represents the a coefficient. 



Instead of assuming a Keplerian accretion flow co-rotating with the central black hole 
as was done in previous work, we include the radial momentum conservation in the equation 
set. I.e., at any radius, the radial velocity is solved by V^T rfl = 0. The angular velocity is 
a little lower than the Keplerian velocity outside about ten gravitational radii to the center. 
In small radii, the radial velocity increases significantly and the angular velocity becomes 
much lower than the Keplerian velocity. The accreted plasma passes through a sonic point 
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at a radius inside the innermost stable circular orbit, becomes super-sonic and finally crosses 
the event horizon of the black hole at a speed of light. We me n tion that we don 't consider 
the causality of the accretion flow (iGammie fc Pophaml Il998l ; iTakahashil 120071 ). We also 
adopt the viscous tensor component derived from thin disk and the a-viscosity coefficient 
for simplicity Such treatments are enough to generate rough dynamical structure for the 
accretion flow, which can be used for observational predictions at 1.3mm wavelength. For 
future observations at much shorter wavelengths and with higher resolution, causal viscosity 
should be calculated carefully to make more realistic images of the vicinity of the black hole. 



r+ ~ 

VIS 



Here, the viscous heating rate is assumed to be equal to the turbulence cascade rate, 
c s B 2 /8irrHg, where B is the magnetic field amplitude determined by the gas pres- 
sure and assumed /3 P , c s is the sound speed, and Hg is the angular half-thickness of the 
disk. The radiative cooling is considered for electrons only and the viscous heating rate for 
electrons is assumed to be different from that for protons, and the radiative cooling for elec- 
trons is much more effective. We adopt th e heating mechanism by turbulent plasma wave 

= a R kBT e T~l, where the 



flBlandford & Eichleri Il987l : iLiu et al. 



20061 ) for electrons as T" 



timescale r acc = 3CirHgv e /cj. 



a, 



x[(3K 3 (x) + K^x)) / AK 2 (x) - 1] with K^x) as the ith 



order modified Bessel function of dimensionless temperature x, and v e is the mean electron 
speed. Therefore the energy equation for electrons becomes 



nu r a P k 



dT^ 
dr 



A 



rad' 



(3) 



and the energy equation for both protons and electrons derived from local energy conservation 



V M T^ = is 



nu r -^[r]S{m p + m e )\ = T+ is - A rad , 



(4) 



where is the Coulomb energy exchange rate between electrons and protons and A~ d is the 
radiative cooling rate calculated from the sum of the energy losses due to the synchrotron 
radiation, the synchrotron self-Comptonization of the soft photons and the bremsstrahlung. 
The temperatures of protons and electrons can be solved with the EqJH and EqJ3] combined. 

The dynamical structure of the accretion flow is determined by solving all the physical 
quantities on the equatorial plane and the scale height from the outermost radius of 10 4 
Schwarzschild radii to the event horizon. We assume homogeneous distribution in the vertical 
direction at any radius and a configuration of magnetic field to be parallel to the velocity 
field in the upper side of the accretion flow and reversed in the lower side, which forms 
from the initial field lines in vertical structure and the shearing of the Keplerian accretion 
flow. Considering synchrotron emission of only thermal population of electrons together with 
self-absorption and birefringence effects, full general relativistic radiative transfer along ray 
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trajectories is then perfor med to give po l arized imag es at any given frequ ency. Corresponding 
equations can be found in iHuang et al.l (120081 ) and iHuang et al.l (120091 ) . 



The model parameters are set in the following proc edure. First of all, (3 U is fixed to 
be 0.7 when the turbulence saturates (IPessah et al.l 120061 ) . Secondly, the inclination angle i 
is fixe d to be 45° as a typical value to r eproduce high linear polarization in sub-millimeter 
band (jLiu et al.ll2007l : IHuang et al.ll2008l ). The position angle is chosen as 140°, which can 
make the accretion flow reproduce the observed position angle of electric vector (EVPA) from 



millimeter to near-infra red band ( iBower et al.l 120051 : iMacquart et al.l 120061 : iMarrone et al. 
20071 ; iMeyer et al.l 120071) with the mean external rotation measure of ~ 5.6 x 10 5 rad • m -2 



adopted (IMarrone et al.l 120071 ) . Next, we arbitrarily set a value for spin a and choose an 
appropriate value for (3 p . In practice, if a is set, [3 P should be small enough to keep very low 
circular polarizations in the sub-millimeter band reported by IMarrone et al.l (120061 ). Finally, 
we adjust values for M and C\ and calculate the synchrotron radiation of ther mal electrons 
to re produce the linear polarization degree of ~ 10% in sub-millimeter band (lAitken et al. 



20001 ) and overall spectrum from centimeter to near-infrared band, fixing the flux density at 
1.3 mm to 2.4 Jy. In practice, the linear polarization is most sensitive to M and the flux 
density is affected by the combination of M, C%, and (3 p . 

We adopt i = 45°, O = 140°, which can explain the polarization observations well, as our 
fiducial model, then choose three sets of model parameters corresponding to three different 
spins. In detail, those parameters are a = 0, [3 P = 0.4, C\ = 0.202, M = 5 x 10 17 g ■ s~ x 
for a non-rotating black hole, a = 0.5, f3 p = 0.2, C\ = 0.357, M = 3 x 10 17 g • s -1 for a 
mildly-rotating black hole, and a = 0.9, (3 P = 0.1, C\ = 0.765, M = 1.2 x 10 17 g • s" 1 for a 
fast-rotating black hole. We show the global solutions of the three cases in Fig{T] and the 
spectra predicted by them in Figj2j The observed emission bump and high linear polarization 
degree in th e sub-millimete r band are predicted by the accretion flow. Another jet/outflow 
component (jLiu et al.ll2007l ). which is not shown in these figures, contributes to the emission 
and depolarization in millimeter and longer wavelengths. For the wavelength of 1.3mm we are 
interested in, we assume the emission is dominated by the accretion flow and the jet/outflow 
component may produce a weak depolarization and a certain change in its EVPA. 



3. Potential (u,v) Coverage For 1.3MM VLBI 

We consider the following eight potential mm- VLBI stations to simulate the 1.3mm 
VLBI observations in this paper: Hawaii (H), including JCMT and Submillimeter Array 
(SMA); SMTO (S); CARMA (C); the Large Millimeter Telescope (LMT, L) on Sierra Negra, 
Mexico; Chilean station (A), consisting of Atacama Submillimeter Telescope Experiment 
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(ASTE) and Atacama Large Millimeter Array (ALMA); the IRAM Plateau de Bure (PdBI, 
P), France; the IRAM Pico Veleta, Spain (PV, V); and a proposed station locating on 
Mount Cook in New Zealand (MT-COOK, M), one of the candidate positions for telescope 
establishment in the future. The (u, v) coverage from these stations is shown in the top-left 
panel of Figj3l These stations offer a (u, v) coverage in almost all the direction and various 
baseline lengths. However, it takes about 11 hours to complete these tracks due to Earth 
rotation. Here we divide the total observing time into six parts, with the interval of ~ 
2hr. In FigJHJ they are shown in square (black), cross(blue), triangle (magenta), circle (violet), 
rhombus (green), and plus(red) symbols, respectively, from the beginning to the end. The 
(u, v) coverage in the same or neighbour symbol, is considered for a simultaneous observation. 

We choose three typical sub-coverages from the total observation. Sub-coverage i is a 
sub-coverage during ~ 16 — 20hr (UT), produced by the Hawaii (JCMT included), SMTO, 
and CARMA, shown in triangle and circle in the top-r ight panel of Fig | 3l Th is sub-coverage 



includes the coverage obtained in the observations by iDoeleman et al.l (120081 ). Sub-coverage 
ii is during ~ 14 — 18hr (UT) from all the stations considered, shown in cross and triangle 
in the bottom-left panel in Figj3j In this sub-coverage, a wide range of baseline length of 
~ — 5.5GA in northwest-southeast direction is included. Especially, lengths of ~ 3 — 5GA 
are covered by baseline group AC/AL/AS. In almost its perpendicular direction, i.e. the 
northeast-southwest direction, comparable lengths are also covered by baselines CH, HS, 
and HL. Simultaneous detections on both perpendicular direction s are very importa nt to 



determine the geometry of the assumed accretion flow, as shown in iHuang et al.l (120071 ) and 
as follows. Sub- coverage in is also from all the eight stations, but during ~ 20 — 23hr (UT), 
shown in rhombus and plus in the bottom-right panel. In this sub-coverage, a wide range of 
baseline length of ~ — 7GA in near east-west direction is covered by baselines CS, CH, HS, 
HL, and AM. The baselines AL, AS, AC, AH, HM, CM, MS, and LM cover a wide range of 
directions apart from the east-west, but only in lengths of > 4GA. 



4. Visibility Analysis And Parameter Estimation 
4.1. Black Hole Spin Estimation For the Fiducial Model 



For the fiducial model described in Secj2j with inclination angle i = 45°, position angle 
= 140°, but with three different black hole spins, a = 0, 0.5, and 0.9, we convolve the 
images with the scattering ellipse. We adopt mean FWHMs of the scattering screen of 
1.39A 2 for the major angular size in milli-arc-second and 0.69A 2 for the minor one, and the 
orient ation of 80° for th e major axis, which is derived from size measurements in centimeter 
band (jShen et al.l 120051 ). We then perform 2-dimensional Fourier transform with special 
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(u, v) coverage provided by those potential stations. With the sampling of sub-coverage i, 
the corresponding visibility profiles are shown in FigJH For each case, the scatter-broadened 
image is shown in the left. Its black hole shadow structure is shown as a white region at the 
center. The dark-red region beside the shadow shows the accretion flow emission boosted 
by Doppler effect. The predi c ted v isibilities are shown in triangl e and circle , with the (u,v) 
coverage in 



Doeleman et al.l (120081 ) included. As mentioned in iFish et al.l (120091 ). Sgr A 



could be considered as observed in its quiescent state since it exhibited the total flux density 
lower than other measurements. These data can be compared with our predictions within 4 
hours, although they were measured on two consecutive days. Their measurements of total 
flux density, correlated flux density on baseline CS, and correlated flux density on baseline 
HS are marked by filled circle, filled square, and filled triangles, respectively. The upper 
limit on baseline CH is marked by a downward arrow. 

The models with a = and a = 0.5 provide good explanations to the visibility measure- 
ments. Notice that the predicted visibilities in triangle at baselines CH and HS, in similar 
directions, cannot be fitted by a Gaussian profile. Clear structures of null point, or valley 
point, are predicted at ~ 3GA, which are due to the existence of the black hole shadow 
(IMiyoshi et al.ll20071 ; iHuang et al.ll20071 ). We note that this structure may be absent in a lot 
of cases (see in the following). In our fiducial model, however, the appropriate lengths and 
directions make these two baselines, CH and HS, potential to detect the important structure, 
if better sensitivity can be achieved for CARMA. 

When the black hole is assumed to be fast rotating, the predicted correlated flux densities 
on baselines CH and HS also increase. This is because that the emission region become more 
compact with large spin. When the black hole spin increases to 0.9, the predicted visibilities 
on baseline CH nearly exceed the observed upper limit. Furthermore, the structure of the 
valley point cannot be recognized. The emission region becomes so compact and totally 
dominated by the blue-shifted side. Therefore, for the fiducial model with % = 45°, G = 140°, 
determined by polarization data, the 1.3mm visibility measurements prefer a black hole spin 
parameter of a < 0.9. 



4.2. Orientation Estimation With 1.3mm VLBI 

The fiducial model with a relatively mildly rotating black hole can reproduce recent 
1.3mm visibility measurements. The orientation of the accretion flow is constrained by 
polarization observations, based on the assumption of a large-scale magnetic field in toroidal 
configuration. However, this specific orientation is not conclusive without direct imaging. In 
spite of the polarizations, the measurements on baselines CH and HS can be explained by 
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many appropriate combinations of a, i, and G, e.g., highly-inclined disk with O ~ 90° — 150° 
for relatively mild spin, or nearly face-on disk with any position angle and spin. Our favorite 
parameters are relatively mildly-rotating black hole, mild inclination angle, and position 
angle toward southeast. On the other hand, any model in opposite parameters, i.e., nearly 
face-on, nearly egde-on, position angle toward northeast, or extremely high spin, can be an 
extreme alternative to our favorite model. 

We focus on the disk orientation, i.e., i and O, first. Choose one fiducial model with 
a = 0.5, % = 45°, O = 140°, and set three alternative cases with a = 0.5, i = 10°, 6 = 140°, 
a = 0.5, % = 90°, 6 = 140°, and a = 0.5, % = 45°, O = 50°. We sample the predicted 
visibilities with the sub-coverage ii mentioned in SecJSJ which covers comparable baseline 
lengths in almost perpendicular directions by baseline groups CH/HL/HS and AC/AL/AS. 
In particular, baseline AL can perform simultaneous observations with baseline HL at ~ 4GA, 
and baselines AL and AS can perform quasi-simultaneous observations (difference within 4 
hours) with baselines CH and HS at ~ 3 — 3.5GA. The corresponding visibilities are predicted 
in Fig|5j shown in cross and triangle. It can be found that the visibilities predicted by the two 
baseline groups make big differences among all the possible cases. We can divide the visibility 
results into three groups: Group I, the predicted visibilities yielded by the two orthogonal 
baseline groups are comparable, e.g., the case shown in the top-right panel; Group II, the 
predicted visibilities yielded by AC/AL/AS are in higher correlated flux density compared to 
those yielded by CH/HL/HS, e.g., cases in the left two panels; Group III, opposite to Group 
II, e.g., the case in the bottom- right panel. These observations are capable of determining 
the basic geometry of the intrinsic emission region of Sgr A*. Detailed discussions are as 
follows. 

If the results in Group I are obtained, the emission region is almost symmetrical to 
the east-west direction. It is probably the nearly face-on case, similar to that shown in the 
top-right panel of Fig0 We show its visibilities by the total (u, v) coverage in the left panel 
of FiglHl In such a case, the visibilities show a similar profile, with a clear valley along 
all directions. Baselines with appropriate length, ~ 2 — 3GA in this case, in any direction, 
such as CH, HS, CL, AL, have the potential to detect such structure. However, there are 
exceptions. Cases with higher inclination angle can be included in Group I, if the position 
angle is close to north or south. We show an example with a = 0.5, i = 45°, O = 165° in 
the right panel of Figj6j In this case, the visibilities again show a similar profile along all 
directions, though for all the baselines we consider here the valley is missing. However, these 
profiles significantly deviate from Gaussian profiles. 

If the results in Group II are obtained, the nearly face-on cases are excluded, and the 
position angle is estimated ranging from ~ 80° to ~ 150°, or from ~ —30° to ~ —100°. 
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The structure of valley point may be detected on baselines CH or HS at ~ 2.5 — 3.5GA, 
except the highly- inclined cases, i.e., i > 60°. In general, for a given inclination angle, the 
differences in correlated flux density between AC/AL/AS and CH/HL/HS are predicted to 
reach maxima at ~ 120°. For a given position angle, such differe nces are predicted to be 



larger with a higher inclination angle. VLBI observations reported in lDoeleman et al.l (120081 ) 
are insufficient to produce a unique estimate of the disk orientation in this case. However, 
Earth-aperture synthesis with additional VLBI stations may be able to do so. Sub-coverage 
Hi is complementary to sub-coverage ii because of the different directions it covers. We 
sample the predicted visibilities with the new sub-coverage for the two cases in left panels 
of Figj5l We then make the same sampling for these two cases by changing their position 
angles into 95°. Predicted visibilities of all these four cases are shown in FigJTl 

Two baseline groups are considered, AM/CH/HL/HS which provides coverage on near 
east-west directions, and AC/AH/AL/AS/CL/LS which provides coverage from northeast 
to northwest. These two baseline groups can perform simultaneous observations with com- 
parable baseline lengths but significantly different directions. In particular, baselines LS and 
CL track ~ 45° from baseline CH ~ 1.5 — 2.5GA. Baseline AL tracks nearly perpendicular 
to baseline HL at ~ 4GA. And AM and AS track nearly perpendicular to baseline AM at 
~ 5 — 5.5GA. With a position angle close to ~ 150° (or ~ —30°), e.g., the two cases with 
6 = 140°, the visibilities yielded by the two baseline groups are predicted to be comparable 
at < 5GA, while in larger difference with higher inclination angle at > 5GA, i.e. between AC, 
AS, HM, and AM. If the position angle decreases from ~ 150° to ~ 80° (or from ~ —30° to 
~ -100°), e.g., the two cases with 6 = 95°, the visibilities yielded by group AM/CH/HL/HS 
are predicted to become much lower than those yielded by group AC/AH/AL/AS/CL/LS. 
With higher inclination angle, these differences become larger, in particular, between AL 
and HL at < 4GA, and between HM (or AS) and AM at < 5 — 5.5GA. In addition, the valley 
structure, or deviation from Gaussian may appear at baselines HL and HS at < 3 — 4GA. 

With the two sub-coverages discussed above, we may make an estimation with an error 
bar as small as ~ 20° of the position angle for an accretion flow of Group II. Furthermore, 
an edge-on case may be easily recognized by visibilities at > 5GA. 

Cases in Group III are not preferred by the recent 1.3mm measurements. Those mea- 
surements have strictly constrained these cases with a < 0.5, 20° < i < 45°, and 9 ~ 0°, so 
that we won't discuss more in this paper. Of course, we can further use sub-coverage Hi for 
confirmation. 
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4.3. Black Hole Spin Estimation With 1.3 mm VLBI 

Small values of black hole spin parameter are preferred by recent measurements. How- 
ever, this is not conclusive because it is based on our fiducial model, with assump t ions o f 



specific electrons heating mechanism and magnetic field structure. In lYuan et all (12 009) 



we ad opted the general radiatively inefficient accretion flow based on the work of iManmoto 



( 120001 ) and I Yuan et all (120031 ). and found the recent measurements prefer extremely high 
values of black hole spin, in stark contrast to the result in this paper. Therefore, the spin 
parameter inferred from observational visibilities depends significantly upon the adopted 
model. 

Given a specific model, the emission region becomes more compact with the increase of 
the black hole spin, so that the structure of the visibility valley, if present, moves towards 
longer baselines. What can be concluded from the visibility measurements is the basic 
geometry, e.g., compactness, or asymmetry, of the emission region. However, determining the 
black hole spin from this is generally model-dependent. We show the images and predicted 
visibilities of two cases for the fiducial model adopted in this paper (M odel A) in the left 



column in FigJHJ For comparison, two cases for the model adopted in I Yuan et all (120091 ) 
(Model B) with the same parameters of spin and orientation are shown in the right column. 
It is found that the emission region of Model A with a non-rotating black hole is even more 
compact than that of Model B with a fast-rotating black hole, e.g., a = 0.9. Notice the 
predicted visibilies on baseline CH and HS, in circle and triangle at < 3 — 3.5GA. and the 
data of recent measurements. Fast-rotating black hole is excluded by Model A, but preferred 
by Model B. 

Such a difference is mainly caused by the different electron heating mechanism used. In 
Model A, the heating mechanism is dominated by turbulence in plasma wave. This heating 
mechanism is so effective that the electron temperature exceeds 10 11 K in the innermost 
region. The angular velocity of the accretion flow is close to the Keplerian velocity, Q > 
0.9fi^, outside the sonic point. The strong Doppler boosting in the blue-shifted side makes 
the emission region more compact. In Model B, the electrons are heated by a small fraction 
of the viscous heating. The electron temperature keeps lower than 10 11 K even near to event 
horizon, except in cases with extremely fast-rotating black hole. Compared to Model A, 
it shows a flat profile in the innermost region, which extends the emission region with the 
same spin. Moreover, the angular velocity of the accretion flow is much lower than the 
Keplerian velocity, Q < 0.6Q^. The much weaker Doppler effect makes the emission region 
more extended. 



In order to make a good estimation of the black hole spin parameter, we should take 
advantage of the properties of the space-time geometry that are independent of configuration 
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of the emission region. For example, the shape, size of the black hole s hadow, and its 
horizontal shift from the rotation axis of the black hole (ITakahashil 12004 ) are the same 
for different models with the same spin parameter, as if the emission region is assumed to 
extend to the event horizon. Therefore, high-fidelity imaging via VLBI at mm/sub-mm 
wavelengths is critical for this. Detecting the innermost stable circular orbit of hot-spot 
structure , if exits, during flaring activity is another way to constrain the spin parameter 
directly faroderick fc Loeblboosthood : boeleman et alJbo09h . 



5. Discussion 

We consider eight potential stations to simulate the future 1.3mm VLBI observations 
for Sgr A* to test the specific accretion flow model with plasma wave heating mechanism. 
These stations provide a highest angular resolution of ~ 30/ias, which is comparable to the 
apparent angular size of the predicted black hole shadow structure. For our fiducial model 
with % ~ 45° and G ~ 140°, the recent 1.3mm VLBI measurements constrain its black 
hole spin parameter as a < 0.9. More detailed visibility analysis suggests that future VLBI 
observations can improve these estimates substantially. 



5.1. Comparisons To Earlier Work 



As mentioned in the previous section, the black hole parameter estimation is dependent 
on the specific model adopted to fit the observations. This is born out by other efforts to fit 
the 1.3mm- VLBI results. 



Broderick et al.l (120091 ) followed the general RIAF model but adopted a Keplerian ve- 
locity distribution. They preferred small values of black hole spin. This is different from 
the prediction of extremely large spin from our Model B, which may be caus ed by different 



treatments of velocity field and gravitational effect, i.e. I Yuan et al.l (120031 ) used Pseudo- 
Newtonian potential instead of Kerr metric. Interestingly, their likeliest values of black hole 
parameter, a = 0, % — 40°, = — 23°(157°), are consistent with those preferred in our Model 
A. The Keplerian, or near Keplerian, velocity field generically makes the emission region 
more compact, preferring small spin parameter. 

Their most probable combination is with values a — 0, i — 90°, = — 18°(162°). This 
orientation is consistent with the preference of a high inclination angle and a small posi- 
tion angle for the s ame RIAF model, constrained by size measurements in 7mm and 3mm 
( iHuang et al.l 120071 ). However, our Model A with the same combination is excluded by ob- 
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servational data because it predicts too compact an image for the emission region. This is 
due to our use of the plasma wave heating mechanism. 

There are other efforts of estimations on the black hole parame ters, especially orien- 



tation, based on different observations and theoretical models. E.g., iMarkoff et al.l (120071 ) 



derived a high inclination angle and a favorite value of ~ 105° for position angle , by compar 



i ng th e visibilities observed in 7mm and predicted by a jet-dominated model; Meyer et al. 



(120071 ) derived a high inclination angle and a range of ~ 60° — 108° for position angle, 
by assuming that an orbiting blob contributes to near-infrared flares and a jet component 
contributes to mm-to-submm emission. These models describe morphologies at 1.3mm wave- 
length totally different from an accretion flow, which can only be distinguished by the future 
direct imaging. 



5.2. Evaluations On The (u,v) Coverage 

In order to obtain a good (u, v) coverage, we need simultanesous observation on baselines 
with various lengths and directions. According to theoretical predictions, baseline lengths of 
< 3GA are very important because the valley point structure of the visibility profile may be 
detected. Fortunately, there are five potential stations, Hawaii (including JCMT), CARMA, 
SMT, LMT, and the Chilean (including ALMA), can provide such suitable baseline lengths. 
Although the coverage is still poor, they can perform simultaneous observations at similar 
lengths but at almost perpendicular directions. The basic geometry of the emission region 
of Sgr A* may be estimated by them. Furthermore, ALMA provides excellent sensitivity on 
baselines as long as ~ 7GA (on baseline AH). It is expected to make a critical contribution 
to future sub- mm VLBI observations of Sgr A*. 

The proposed station in New Zealand, MT-COOK, also makes great contribution. It 
provides long baselines at a large range of direction. With its participation, these eight 
stations cover almost all the directions at lengths of > 5GA. The longest baseline length 
of ~ 8GA, which can provide the highest resolution, is also achieved by MT-COOK (on 
baselines CM and MS). 

With the Earth's rotation, these eight stations produce sufficient (u, v) coverage to 
produce images of Sgr A*. However, it takes almost 11 hours to complete. In this paper, 
we treat observations within 2 hours as the simultaneous and those within 4 hours as the 
quasi-simultaneous. The estimate of the basic geometry of the emission region is made from 
the hypothetic observation of four hours duration, e.g., with sub-coverage i mentioned in 
SecfS] adopted. Thus, we have to assume that Sgr A* is in its stable state during the entire 
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4-hour period. If we want to make more detailed estimates using coverage in a larger range 
of directions, we must assume that this stable state exists even longer. However, a number 
of flares have been observed in Sgr A* at 1.3mm, on timescales as short as 30 minutes. 
Therefore, reliable estimates can only be made when either Sgr A* is quiescent, the flaring 
region is subdominant at 1.3mm, or the observation lasts for < 1 hour (resulting in poorer 
coverage and less opportunity to detect the visibility valley associated with the black hole 
shadow) . 

There is one more short-coming in the (u, v) coverage. Region with baseline lengths of 
> 5GA is covered in almost all the directions. However, the most important region with 
baseline lengths of > 3GA is lack of coverage in northeast-southwest direction. During any 
simultaneous observation within 2 hours, coverage in two perpendicular directions can be 
obtained. However, at least one more direction between those two is necessary for a good 
estimation of the emission region geometry, also for a higher possibility of detection on the 
valley point structure. This can only be improved with an additional new station included, 
e.g., one in New Zealand. 
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Fig. 1. — Global solutions of the accretion flow with plasma wave heating mechanism. The 
solid lines represent the results for a = 0,/3 p = 0.4, Ci = 0.202, M = 5 x 10 17 g • s -1 . The 
dash-three-dotted lines are for a = 0.5, (3 P = 0.2, C x = 0.357, M = 3 x 10 17 g • s -1 . The dotted 
lines are for a = 0.9, (3 P = 0.1, Ci = 0.765, M — 1.2 x 10 17 g • s -1 . Top-left: radial velocity 
and sound speed; Top-middle: angular momentum, with the Keplerian one in solid grey line; 
Top-right: ratio of scale height to radius; Bottom-left: temperatures of protons and electrons; 
Bottom-middle: number density of electrons; Bottom-right: amplitude of magnetic field. 
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Fig. 2. — Spectrum of three accretion flow models with dynamical structures shown in FigJTJ 
i = 45°, and = 140°. Top: spectrum of total flux density of synchrotron emission; Middle: 
linear polarization degree of synchrotron emission; Bottom: position angle of electric vector 
(EVPA) with an external rotation measure of ~ 5.6 x 10 5 rad • m -2 considered. 
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Fig. 3. — Top-left: (u, v) coverage produced by eight potential stations at 1.3mm wavelength. 
H: Hawaii, including JCMT and SMA, S: SMTO, C: CARMA, L: LMT, A: Chilean, including 
ASTE and ALMA, P: PdBI, V: PV, M: MT-COOK. Tracks in different parts of observational 
time are shown in different symbols (square, cross, triangle, circle, rhombus, and plus), with 
the interval of ~ 2hr. Top-right, sub-coverage % during ~ 16 — 20hr (UT), produced by H, 
S, and C. Bottom-left: sub-coverage ii during ~ 14 — 18hr (UT), produced by all the eight 
stations. Bottom-right: sub-coverage Hi during ~ 20 — 23hr (UT), produced by all the eight 
stations. 
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Fig. 4. — Images and visibilities predicted by the fiducial model with a = 0, 0.5, and 0.9 , 
with the sampling of sub-coverage i. The observational data are from lDoeleman et al.l (120081 ). 
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Fig. 5. — Images and visibilities predicted by four cases with a = 0.5 but different orien- 
tation, with the sampling of sub-coverage ii. In these four cases, the visibilities yielded by 
baseline group AC/AL/AS can be higher than, lower than, or comparable to the visibilities 
yielded by baseline group CH/HL/HS, see in context for details. 
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Fig. 6. — Images and visibilities predicted by two cases with a = 0.5 but different orientation 
combinations, with the sampling of the total (u, v) coverage. In both the cases, the visibilities 
yielded by baseline group AC/AL/AS are comparable to the visibilities yielded by baseline 
group CH/HL/HS. 
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Fig. 7. — Images and visibilities predicted by four cases with a = 0.5 but different orientation 
combinations, with the sampling of sub-coverage in. In these four cases, the visibilities 
yielded by baseline group AM/CH/HL/HS can be higher than, lower than, or comparable to 
the visibilities yielded by baseline group AC/AH/AL/AS/CL/LS, see in context for details. 
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Fig. 8.— Left: Images and visibilities predicted by the fiducial model adopted in this 
paper (Model A) with a = and 0.9, with the sampling of the total (u,v) coverage. Right: 
Images and visibi lities predicted by radiatively inefficient accretion flow model adopted in 
Yuan et al.l (120091 ) (Model B), with the same spin, orientation, and (u,v) sampling as the 
left ones. 



